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We present a noneqmlibrium nonperturbative field theory for the Kondo effect in strongly in- 
teracting quantum dots at finite temperatures. Unifying the slave-boson representation with the 
Keldysh field integral an effective Keldysh action is derived and explored in the vicinity of the zero 
slave-bosonic field configuration. The theory properly reflects the essential features of the Kondo 
physics and at the same time significantly simplifies a field-theoretic treatment of the phenomenon, 
avoiding complicated saddle point analysis or 1/N expansions, used so far. Importantly, our theory 
admits a closed analytical solution which explains the mechanism of the Kondo effect in terms of 
an interplay between the real and imaginary parts of the slave-bosonic self-energy. It thus provides 
a convenient nonperturbative building block, playing the role of a "free propagator", for more ad- 
vanced theories. We finally demonstrate that already this simplest possible field theory is able to 
correctly reproduce experimental data on the Kondo peak observed in the differential conductance, 
correctly predicts the Kondo temperature and, within its applicability range, has the same universal 
temperature dependence of the conductance as the one obtained in numerical renormalization group 
calculations. 



PACS numbers: 72.15.Qm, 73.63.-b, 72.10.Fk 
I. INTRODUCTION 

Experimentally discovered almost eighty years agcP 
and qualitatively explained thirty years after^ , the Kondo 
effect, a minimum in the temperature dependence of the 
resistance in magnetic alloys, was later anew brought 
into the world and garbed with physics of quantum dots 
(QD)'^'*. That time this complicated phenomenon was 
supplemented by physics of nonequilibrium due to the 
dot coupling to contacts biased by an external voltage. 
In this setup the Kondo effect, a nonperturbative phe- 
nomenon induced by both the electron-electron interac- 
tions and the QD-contacts coupling, appears as a sharp 
many-particle resonance in the tunneling density of states 
(TDOS) at the Fermi energy. An immediate consequence 
of this resonance is the zero-bias maximum observed ex- 
perimentally in the QD differential conductance at low 
temperatures. 

Theoretical prediction^^lEl of this behavior made be- 
fore the actual experiments use the single-impurity 
Anderson model (SIAM^P with the on-dot interaction 
U. Those early works utilized either a quasi-particle 
transformatiorP to analytically predict resonant trans- 
mission through QDs for arbitrary U or the noncross- 
ing approximati on (N C A) , equations of motion and per- 
turbation theorjlSEHHm to numerically and analytically 
describe several important stationary and nonstationary 
features of the Kondo effect in strongly interacting QDs 
modeled by the infinite-C/ SIAM, in particular, in a slave- 
boson representatiorPsHHI 

Later, as an alternative to the slave-boson representa- 
tion, diagrammatic expansions within the reduced den- 
sity matrix formalisrrP^ were applied to SIAM for the 
infinite-?/ case. Extracting a certain infinite subset of 
diagrams, an analytical nonperturbative expression for 
the TDOS was obtained at finite temperatures. Here the 
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FIG. 1: (Color online) Classes of the slave-bosonic theo- 
ries for the Kondo effect in QDs. Their range of quantita- 
tive reliability with respect to the Kondo temperature Tk 
(Tk/F ~ exp[-27r(/io - ed)/F], where F is the total QD- 
contacts coupling strength, which is twice that of Ref. fio 
the QD chemical potential, ed the single-particle energy level 
of the QD) is given by their horizontal location. The verti- 
cal location of a class shows its type, numerical, analytical or 
both. As to our field theory, what the blue class shows is what 
one generally expects from nonperturbative field theories in 
the vicinity of the zero slave-bosonic field configuration. 



absence of the double occupancy played a crucial role in 
identifying the relevant diagrams. 

For small, intermediate and large U the Kondo prob- 
lem in QDs was widely explored numericalljff^Hini 

In the present paper we develop an analytical nonequi- 
librium real-time field theory of the Kondo effect in 
strongly interacting QDs at finite temperatures using the 
infinite-f7 SIAM. The field-theoretic approach is based 
on the slave-boson representation and the Keldysh field 
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integraPil^. 

Importantly, our theory is nonperturbative in both the 
electron-electron interaction and QD-contacts coupling. 
The necessity in such a theory is obvious from its advan- 
tages. Firstly, it is a field theory and, thus, it has a clear 
systematic generalization from the present relatively sim- 
ple basic model to more involved setups such as the ones 
with ferromagnetic contacts, superconducting contacts, 
finite U systems, etc. This is especially important since 
there is a limited number of analytical theories nonper- 
turbative in both the electron-electron interaction and 
QD-contacts coupling and which at the same time have 
a straightforward generalization scheme. Secondly, it is 
an analytical theory and, thus, it will help to reveal a rel- 
evant physical picture behind new and more complicated 
physical systems as, e.g., those mentioned above. This 
is definitely an advantage over numerical methods such 
as, e.g., NCA, which provide good quantitative descrip- 
tion but hide the essence of the physics, giving only an 
indirect access to it. 

In Fig. [ijwe show a comparative layout of the slave- 
bosonic theories for the Kondo effect in QDs. A basic 
nonperturbative analytical slave-bosonic field theory of 
the Kondo effect in QDs within the range from below 
Tk to higher temperatures, i.e., within the most relevant 
experimental temperature range (blue (uppermost) class 
in Fig. [T]), is formulated in this work. 

So far, using the slave-boson method in the con- 
text of QDs, th e resul ts were obtained nonperturba- 
tively numericall} *^ ! '' 1 ^ 1 ^ ^ below and/or above Tk, pertur- 
batively semi-analyticall}ff^ above Tk and nonperturba- 
tively analytically2ll below Tk . 

The slave-boson approach excludes the double occu- 
pancy of a QD restricting its Hilbert space to zero and 
single occupancy. This restrictio n was taken into account 
exactly in the numerical solution j^*^^^ * using an additional 
integral removing the constraintP^. However, since the 
mean fielcpSl and 1/N expansiorPS theories, using the 
same integral trick, impose the constraint only approxi- 
mately, they fail at higher temper at ured^. 

Here our goal is the development of an analytical 
ground for the slave-bosonic nonperturbative field the- 
ories which could be placed within the blue (uppermost) 
class in Fig. [T] To avoid problems with high temper- 
atures we do not use the integral trick to take into ac- 
count the restriction of the Hilbert space but instead use 
an alternative method based on taking a certain limit 
with respect to a real parameteil^Sl. This alternative 
method has not been applied so far in conjunction with 
the Keldysh field integral. We demonstrate that the cor- 
responding limit can be exactly performed analytically 
after the Keldysh field integral for the TDOS has been 
analytically calculated using a certain approximation. 

Specifically, the approximation concerns the effective 
Keldysh action obtained after integrating out all elec- 
tronic degrees of freedom. This action being a nonlinear 
functional of the slave-bosonic field is expanded up to 
second order in this field. Note, that this does not im- 



ply any perturbation since the action itself is the argu- 
ment of an exponent. According to the general concept 
of the condensed matter field theorjfSH the physics of our 
model is the physics in the vicinity of the zero slave- 
bosonic field configuration. We demonstrate that this 
physics contains the Kondo effect in QDs at finite tem- 
peratures. This scenario complements^- the saddle point 
analysia^l- valid deep below the mean field theory slave- 
bosonic phase transition, i.e., at low temperatures when 
the slave-bosonic field in the mean field theory is con- 
densed. Surprisingly, our simple or "bare" theory, which 
might play a role of a nonperturbative "free propagator" 
for more advanced nonperturbative theories of the blue 
class in Fig. [T] already provides a good description of the 
Kondo peak observecF' in the differential conductance at 
temperatures close to Tk- 

The paper is organized as follows. In Section [IT] we for- 
mulate the problem on the operator level. Then in Sec- 
tion lllll we translate this formulation into a field-theoretic 
framework using the Keldysh field integral and derive an 
analytic expression for the TDOS. Finally, we discuss our 
results and make conclusions in Sections IV and |Vj re- 
spectively. 



II. OPERATOR FORMULATION: 
HAMILTONIAN AND OBSERVABLES 

We start with the infinite-fj Anderson Hamiltonian in 
a slave-boson representation. As is well known^^, in the 
case when U — oo, the Anderson Hamiltonian, H^i = 
J2a- '^dnd,a + Uhd,^nd^i, where nd,a = dlda, a =t, i, takes 
the form = X^o- ^d/l/cr, where the new fermionic op- 
erators are related to the original ones as d^- — fab\ 
d\ = f\b, and h, are the annihilation and creation 
operators of a slave-boson. Using these new fermionic 
and slave-bosonic operators the tunneling Hamiltonian, 
describing the coupling of the infinite-?/ Anderson QD to 
contacts, can be written as 



(1) 



where Ca, cjj are the annihilation and creation operators 
of the contacts fermions, a includes the contacts complete 
set of quantum numbers and the contacts labels, left (L) 
or right (R), and Taa are the tunneling matrix elements. 
The contacts are described by the Hamiltonian He = 
^aC^Ca and are assumed to be in equilibrium with the 
chemical potentials /xl,r = ^q — bVl.r, with F = Vl — 
being an external voltage. The total Hamiltonian \s H = 

The restriction of the QD Hilbert space to zero and 
single occupancy requires the total number of the new 
fermions and slave-bosons, Q = b^b + Y^^flfa, to be 
equal to one, Q — I. This restriction must be taken into 
account in a QD observable (O). There are two ways 
of doing thalpS. In the QD context only the first one. 



3 



i.e., the integral way was used so fajEEMSlMl^ Here we 
employ the second way from Ref. [25l 



lim e 

/I— )-oo 



lim ef^^'^I[poQe'l^^'Q] 



(2) 



where C/t,t', is the evolution operator with respect to 
the Hamiltonian H, po — exp[— /3(^fd — A*o^d)] (X" 
exp[—f3{Hc^J2x f^xNx)], {x = L,R) is the initial statisti- 
cal operator with and being the number operators 
of the QD and contacts and /3 is the inverse temperature. 
Eq. (H may be rewritten in two equivalent forms, 



and 



lim e'^^Tr[C/_oo,ooC^oo.tOf/t,-ooPoe 
lim e'3MTr[poQe-'3'^Q] 



(3) 



lim e^f'Tr[U^o.,tdUt,ooU^. 

lim e/5''Tr[poQe-^'3A'Q] 



(4) 



These two forms have different interpretation: in Eq. ([s]) 
the observable is taken within the evolution from — oo to 
oo while in Eq. Q it is taken within the evolution from 
oo to — oo. 

To develop a Keldysh field theory one first equivalently 
rewrites Eqs. ^ and Q on the Keldysh contour Ck- 
To this end one gives the creation and annihilation op- 
erators a formal temporal argument to allow the time- 
ordering operator to appropriately interlace operators. 
Afterwards one may take the half-sum of the two equiv- 
alent expressions to get the following symmetric form: 



(o>(t) 



lim 



rTr [Tc^e •'^k 



7VoM->ooTr[/5;)(/x)] 

d(i+) + 0(L)„ 



(5) 



where t+ and t_ are the projections of t onto the forward 
and backward branches of Ck , 



1 _ Tr(pc) 



H' = H + nQ, 



(6) 



lim {eMP^^)^AQpoi^^)]} 

/I— f oo 



and PC is the statistical operator of the contacts. The 
expression under the limit in Eq. ([5| can be written 
as the Keldysh field integrat^l'2^ for a fixed value of p,. 
The basic steps in the construction of the Keldysh field 
integral are identical to the ones presented in Ref. [211 
The details specific to our application of this field integral 
are given in the next section. 

We would like to emphasize the absence in Eq. ([s]) 
of any prefactor depending on the QD-contacts coupling. 



This is a great advantage of the Keldysh field theory 
over non-field-theoretic approaches- and imaginary-time 
field theorie^^H Indeed, in our approach there is no 
need for an independent calculation of such prefactors. 
This greatly simplifies the analytical exact projection 
onto the physical subspace. This fact was not realized in 
the first attemplPf to combine the Keldysh field theory 
and slave-boson approach and, as a result, this attempt 
was reduced to a nonequilibrium analog of equilibrium 
imaginary-time field theories with only an approximate 
projection onto the physical subspace. 



III. FIELD-THEORETIC SOLUTION: KELDYSH 
FIELD INTEGRAL 

We obtain the effective Keldysh field theory in a way 
similar to the one used in Refs. I22|26l for Coulomb- 
blockaded QDs. Namely, we first integrate out the QD 
and contacts Grassmann fields. After this step the field- 
theoretic description is given in terms of the effective 
Keldysh action, S,s[x'\t),xHt)] = S^°\x'\t) , xHt)] + 
Stun[x'^^{t)iX'^{t)]j where x'^''''(^) are the classical and 
quantum components of the slave-bosonic complex field, 
which is just a bosonic coherent stat^^l^ [x'^'(*); x''(*)] 
is the free slave-bosonic action with the standard matrix 
form in the Keldysh space and Stun[x^^{t)iX'^{t)] is the 
slave-bosonic tunneling action, 

Stun[x'Ht),x''m = ~ihtTln[l + TG<^"^, (7) 

where the trace and matrix product are taken with re- 
spect to the temporal arguments and both single-particle 
and Keldysh indices. In Eq. ([t]) the matrices G'-"^ and T 
have the block form in the dot-contacts space. 







G^c\a,t\a't') 





MT{at\a't') 



Ml{at\a't') 




(8) 



(9) 



where the blocks G'f l,{at\a't') are the standard fermionic 
Keldysh Green's function matrices (2x2 matrices in the 
Keldysh space) of the free QD {a = cr) and contacts 
(a = a), 



(0)K. 



GS^-(ai|«'i'), 



(10) 



with G|j^^^' '^{at\a't') being the retarded, advanced 
and Keldysh componentJ^, respectively, and the block 
MT{at\<Tt') is the tunneling matrix (2x2 matrix in the 
Keldysh space), 

M.(atK)4^(^-t')T..i^g:{;J p§^. (11) 
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Note that the only formal difference of Eq. (Ill from the 



corresponding expression in Ref. |22j is that here instead 
of the bosonic phase field we have the slave-bosonic field 
and that this field is not exponentiated. 

Our goal is to investigate a strongly interacting QD 
in the Kondo regime when the Kondo resonance is not 
strong yet. This is the case, e.g., at temperatures above 
the Kondo temperature when the Kondo resonance is 
already present but not fully developed. These range 
of temperatures is most relevant for both experiments 
and practical applications in devices. As the formation 
of the Kondo resonance is related to the QD population 
oscillations induced by the QD-contacts coupling, the QD 
empty state will weakly fluctuate and its probability will 
be small in the regime which we are interested in. Since 
the QD empty state is described by the slave-bosonic 
complex field, those weak oscillations, induced by the 
QD-contacts tunneling coupling, can be obtained from 
the expansion of the tunneling action around the zero 
slave-bosonic field configuration. This expansion is valid 
for small amplitudes of the slave-bosonic complex field 
and thus for small probabilities of the QD empty state. 

We then expand the tunneling action ([t]) around the 
zero slave-bosonic field configuration, i.e., we keep in 
this action only the first non-vanishing term which is 
quadratic in the slave-bosonic field. Additionally, we as- 
sume a Lorentzian contacts density of states, vc{^) — 
vcD{e), D{e) = W^/{e^ + W^), a = {x,k,a} and 
Txkaa' = Sa(T'T. The exprcssion for the tunneling ma- 
trix elements, in particular, means that we, for simplicity, 
consider the case of a symmetric coupling to the contacts. 
The tunneling action then becomes 



c ci 



[x'\t),x''it)] = {g/h)jdtjdt' {r\t) xm) X 

s-(t-i')Wx^'(t')' ^^^^ 



In Eq. (12 \ g = 27r^i^c|Tp and the retarded, advanced 
and Keldysh slave-bosonic self-energies are 



E±(t - t') ^ (z/2) Y^igfit' - t)g^{t - t')^ 

X 

+ g^{t'-t)gfit^t% 

EK(t - t') ^ (^/2) J2{gf{t' t)gf{t t'y 



(13) 



- [aiit t') g-{t - t%g+{t' -t)- g-{t' ~ t)]}, 

where the retarded, advanced and Keldysh components 
of the Green's function matrix are 



de h 

2^ ^^'hu-e±iO' 



'K(hhJ — 6(1 — /i ± zO) ■ 



9x {^) = [g {^) - 9 (^)] tanh 



I3{hu} - fix) 
2 



(14) 



(15) 



with nd(e) being the QD distribution of the new fermions. 

The QD TDOS is defined through the imaginary 
part of the QD retarded Green's function, i^ai*^) = 
-il/hn)Im[G+^^ie)]. 

The Keldysh field integral expression for i'cr(e) is 

X / 2?[x'=\t),x''(0]e*'^''"'''"^*^'''''^*^'x 



dteT-^W 



(16) 



X [x-(i)x+(0)G>(f7t|aO) - x+(i)x-(0)G<(ai|aO)], 



where x± {i) is the slave-bosonic field on the forward and 
backward branches of Ck, 



iG>'<((7t|(T0) = 



11. 

-±--nd(ed + /^) 



exp 



(17) 



As in Refs. I22I26I the distribution nd(e) is the double 
step. In the present context this is due to the noninter- 
acting nature of the new fermions. 

Since the effective Keldysh action Scs[x'^^it),X^{t)] is 
quadratic, the functional integral in Eq. ( 16 ) may be 
performed exactly for any real /i > 0. The limit /i — oo 
is readily taken afterwards. The final nonperturbative 
result for the QD TDOS is 



V{e) 



[ed-e+(.9A)ER(e)]2 + [(g/S)Si(e)P' 



(18) 



where ^^(e), Ei(e) are the real and imaginary parts of 
the retarded slave-bosonic self-energy for which we find 
the following analytical expressions. 



I]R,(e) = h[D{e)/^]{€/2W+ 
+ {l/2T:)KcY[{l+ie/W)tl^[l/2-ip{iW+ fix) 

X 

-^[l/2-i^(^,-e)/27r]]}, 
^i{e)^h[D{e)/2^]Y^nx{e), 



(19) 



where ip{x) is the digamma function and nx{€) are the 
contacts Fermi-Dirac distributions, and the numerator is 



V{e) = 



gD{e) nL(ed) + '^R(ed) ~ 2nL(ed)»^R(ed) 



27r2 



1 - "L(ed)«R(ed) 



(20) 



The imaginary part in Eq. ( |19[ ) has a particularly clear 
physical meaning. The slave-boson describes the empty 
state of the QD. From Eq. ( 19 1 it follows that for energies 



above the chemical potentials fix the imaginary part of 
the slave-bosonic self-energy Si(e) — >• which means that 
the life-time of the empty state goes to infinity, i.e., the 
QD is empty. In contrast, below fix the imaginary part 
is finite leading to a finite life-time of the empty state, 
i.e., the QD is filled. 
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It is important to emphasize that the QD TDOS, Eq. 



(181, is finite at any finite temperature but logarithmi- 



cally diverges for T = at the chemical potentials. Thus, 
one expects that for very low temperatures the simple 
quadratic theory is not valid. This is in accordance 
with what one generally expects on purely theoretical 
ground^^ from expansions around zero field configura- 
tions. Zero field expansions complement low temperature 
theories, being expansions aro und n on-zero field configu- 
rations, like mean field theoriePlal. These non-zero field 
expansions fail to describe the high-temperature behav- 
ior of the Kondo effect in QDs as we have demonstrated 
in Fig. [T] Therefore, one concludes that for a compre- 
hensive description of the Kondo effect in the whole tem- 
perature range it is desirable to have in one's disposal 
both types of field theories. 

Let us clarify the nature of the approximation used 



for the tunneling action, Eq. (12), and the applicability 
range of the resuh for the TDOS, Eq. ((Tsl). Formally 



the small parameter of the expansion is the tunneling 
matrix element T. Hence, F = 2(;/7r should be small. 
However, the tunneling matrix elements always enter in 
product with the slave-bosonic fields. Thus our theory 
will be valid in physical situations where the large values 
of the slave-bosonic amplitude are not important, i.e., 
when the probability of the QD empty state is not large. 
On one side this happens when /io — Cd is large and on the 
other side when the temperature is not too low and thus 
the Kondo resonance is not too strong so that the fluctua- 
tions of the QD empty state are not too large. Therefore, 
we assume that our theory should be applicable for 



/^o — fd ?J r and T >Tk- 



(21) 



One should keep in mind that this is a crude estimate 
and our theory may still work qualitatively and per- 
haps semiquantitatively outside those inequalities. As 
demonstrated in the next section, our simple quadratic 
approximation, accounting for the physics in the vicin- 
ity of the zero slave-bosonic field configuration, provides 
a reasonable description of the Kondo physics. It also 
gives a finite single-particle resonance. This situation is 
definitely better than the one taking place in perturba- 
tive approacheJi^ which have divergences at the single- 
particle resonance. Our theory being nonperturbative 
avoids this problem but requires the next order term, i.e., 
the one which is quartic in the slave-bosonic field, in the 
expansion of the tunneling action ([?]) for a better quan- 
titative description. This will be the subject of a more 
advanced theory. Here we only note that this advanced 
theory may be constructed in terms of the quadratic the- 
ory presented in this work. Namely, the present theory 
will play a role of a nonperturbative "free propagator" 
for the quartic theory. 

It is interesting to look at the result for the QD TDOS, 
Eq. (18), from the physical point of view. Within the 



range of its applicability it suggests that the quasiparticle 
state in the QD represents a superposition of the bare 
electron state and the empty state of the QD coupled 
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FIG. 2: (Color online) The equilibrium and nonequilibrium 
TDOS in the Kondo regime. Here kT = O.OOSF, fio - <^d = 
1.876r, W = lOOr. In equihbrium, V = 0, there is a sharp 
many-particle peak at the Fermi energy. In nonequilibrium, 
V 0, the peak is reduced and split into two lower peaks 
(see inset) . The total spectral weights of the equilibrium and 
nonequilibrium situations are almost the same differing by 
1% which is due to the accuracy of the numerical integration 
involved in the total spectral weight. 



to contacts. Thus in this range the life-time of the QD 
quasiparticles can be estimated as the life-time of the 
QD empty state. This may be very attractive for the 
experiments measuring the quasiparticle life-times since 
observing QD states is easier than a direct observation 
of its quasiparticles. 

Finally, we would like to emphasize an additional fun- 
damental advantage of our method: our field-theoretic 
approach is a tru ly field theory in contrast with mean 
field theorie423124| r^^ye point is that here we work only 
with a physical fleld, being the slave-bosonic field, and 
not with artificial fields like the Lagrange multipliers 
fields used in mean field theories. Thus our theory has 
more transparent access to physics avoiding such arte- 
facts of mean field theories as slave-bosonic condensation. 



IV. DISCUSSION OF THE RESULTS 



In Fig. [2] the QD TDOS ^ is shown. In addi- 
tion to the single-particle resonance at the renormalized 
QD noninteracting energy level it reveals a many-particle 
peak at the Fermi energy, known as the Kondo resonance. 
The formation of the Kondo resonance is explained in 
Fig. [3] as an interplay between SR(e) and Si(e) from 
( [Tq] ). The distance between e — Cd and SR(e) reaches a 
minimum where ^^(e) has its maximum. This happens 
at the Fermi energy. At the same time Si(e) has a steep 
decrease at the Fermi energy. Therefore, the two terms 
in the denominator of Eq. (18) become both minimal at 



the Fermi energy giving rise to the Kondo resonance. 
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FIG. 3: (Color online) The mechanism of the Kondo reso- 
nance formation in the nonperturbative Keldysh field theory 
in the vicinity of the zero slave-bosonic field configuration. 
Here V = Q,kT = O.OOSF, - ed = 1.95r. 
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FIG. 4: (Color online) The Kondo peak in the differential 
conductance. Inset shows the temperature dependence of the 
differential conductance maximum at V = 0. The solid line 
is the result of our nonperturbative Keldysh field theory. The 
circles show the experimental data of Ref. [3] 
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FIG. 5: (Color online) Comparison of our field theory 
with the numerical renormalization group theory. The solid 
line is obtained from the numerical renormalization group 
calculationPJf^. The circles show the differential conduc- 
tance maximum obtained from our field theory for the values 
of the parameters used in Fig. |4] 
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FIG. 6: (Color online) Keldysh field-theoretic prediction of 
the Kondo temperature dependence on the QD single-particle 
energy level (/io - CcO/F. 



To verify our field-theoretic description of the Kondo 
physics we first compare it with experimental data. The 
presence of the Kondo resonance in the QD TDOS has 
an impact on the other QD observables. In particular, 
experiment J^I show a peak in the differential conductance 
at y = 0. Using the expression for the current (see Eq. 
(3) in Ref. ED through a QD together with Eqs. (fl8|) and 



(191 we obtain this behavior of the differential conduc- 
tance shown in Fig. [4j To get Fig. [4] we have taken the 
values of the parameters, F = 2.6875 meV, = 5 eV, 
T = 50 mK, F/[7r(^o - ed)] = 0.1224, close to the ones 
which were estimated in Ref. i3, 

We would like further to compare our Keldysh field 
theory with existing theoretical approaches, in partic- 
ular, with the numerical renormalization group theory 
from Ref. [27] wh ich w as successfully employed to de- 
scribe experiment^^Hlll on the Kondo effect in QDs. As 
we have argued that our theory must be valid for temper- 



atures T ^ 7k , within this temperature range the differ- 
ential conductance maximum in our theory must have the 
same temperature dependence as the one in Refs. 28 2^1 
(see, e.g., Eq. (2) in Ref. [28l where we take s = 0.21) 
for the case of a symmetric coupling. From this high- 
temperature comparison, for each value of /ig — ed used 
to calculate the differential conductance maximum in the 
Keldysh field theory, we can fix the value of Tk to be used 
in the empirical form of Ref. [55] For example, for the 
parameters used in Fig. [4]we get the temperature depen- 
dence of the differential conductance maximum shown in 
Fig. [5| The Kondo temperature Tk = 11 mK agrees with 
the rough estimate Tk < 50 mK given in Ref. [S] where 
a weak asymmetry in the capacitance of the QD to the 
left and right contacts was assumed. The same compar- 
ison between the Keldysh field theory and the numerical 
renormalization group theory can be done for any value 
of /io — Cd for which the Keldysh field theory is applica- 
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FIG. 7: (Color online) Comparison of the universal tempera- 
ture dependence of the differential conductance maximum in 
our field theory and in the numerical renormalization group 
theory. The solid line is the result of the numerical renormal- 
ization group calculationJ^^f^. The circles show the result 
obtained from our field theory. 



ble. In this way we determine how in our Keldysh field 
theory Tk depends on ^LQ — td- This dependence is shown 
in Fig. [6] and it is in full accordance with the standard 
expressioiil^, 



r 



exp 



-27r 



r 



(22) 



where our definition of F is twice that of Ref. [6] (see also 
the caption of Fig. [T]). This proves that our theory for 
the Kondo effect in QDs correctly predicts the Kondo 
temperature. Moreover, our theory, within its applica- 
bility range, also predicts that the differential conduc- 
tance maximum has a universal temperature dependence 
with the scaling given by the Kondo temperature, Eq. 
( p2| ) . This universal temperature dependence is shown in 
Fig. [7]and additionally proves that our Keldysh field the- 
ory, within its applicability range, correctly describes the 
Kondo physics in QDs. As one can see from Fig. [7] the 
Keldysh field-theoretic description is quantitatively reli- 
able for temperatures T ^ 2T'k, which perfectly agrees 
with our theoretical prediction made above in Section|TlT| 



Finally, we would like to say a few words about the 
numerical consistency of our theory. To do this, we em- 
ploy the sum rule given by Eq. (39) of Ref. El In NCA 



this sum rule is always satisfied within 0.5%. In our the- 
ory this depends on how well the applicability criteria, 
Eq. (21), of the Keldysh field theory are satisfied. For 



example, for the parameters presented in Fig. |2]the sum 
rule is satisfied within 15% while for /iq — = 2.5r it 
is 8.5% and for - = 4. OF it is 4.4%. One should 
note that the sum rule is an integral estimate over the 
whole energy range. Thus, the error is gained over the 
whole range of energies. At the same time for a given en- 
ergy the QD TDOS may have higher accuracy. Since the 
only approximation was the truncation of all the terms 
of higher orders than the terms quadratic in the slave- 
bosonic field, to improve the consistency of the method 
and extend its applicability criteria one should go beyond 
the quadratic approximation and this will be done in a 
subsequent study. 



V. CONCLUSION 

We have developed a basic slave-boson nonperturba- 
tive Keldysh field theory for the Kondo effect in quantum 
dots. The theory deals with the physics in the vicinity 
of the zero slave-bosonic field configuration where, as we 
have shown, the main fraction of the Kondo physics is 
located at experimentally relevant temperatures. The 
presented theory has a closed analytical solution for the 
quantum dot tunneling density of states and, despite 
being relatively simple, properly describes experimen- 
tal data on the Kondo peak observed in the differential 
conductance, correctly predicts the Kondo temperature 
and, within its applicability range, has the same univer- 
sal temperature dependence of the conductance as the 
one obtained in numerical renormalization group calcu- 
lations. Therefore, it represents a convenient basis, as 
a free nonperturbative propagator, for more advanced 
theories which could extend the applicability of our ap- 
proach to larger values of the slave-bosonic amplitude 
and, thus, to temperatures much lower than the Kondo 
temperature. 
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